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ABSTRACT 

The sizes of interstellar grains are widely distributed, ranging from a few 
angstroms to a few micrometers. The ultraviolet (UV) and optical extinction 
constrains the dust in the size range of a couple hundredth micrometers to sev¬ 
eral submicrometers. The near and mid infrared (IR) emission constrains the 
nanometer-sized grains and angstrom-sized very large molecules. However, the 
quantity and size distribution of micrometer-sized grains remain unknown as they 
are gray in the UV/optical extinction and they are too cold and emit too little 
in the IR to be detected by IRAS, Spitzer, or Herschel. In this work, we employ 
the ~ 3-8 yum mid-IR extinction which is flat in both diffuse and dense regions to 
constrain the quantity, size, and composition of the /im-sized grain component. 

We hnd that, together with nano- and submicron-sized silicate and graphite (as 
well as PAHs), pm-sized graphite grains with C/HRil37ppm and a mean size 
of ~ 1.2 pm closely fit the observed interstellar extinction of the Galactic diffuse 
interstellar medium from the far-UV to the mid-IR as well as the near-IR to mil¬ 
limeter thermal emission obtained by GORA'/DIRBE, GGRR/FIRAS, and Planck 
up to A < 1000 pm. The pm-sized graphite component accounts for ~ 14.6% of 
the total dust mass and ~2.5% of the total IR emission. 
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1. Introduction 


The absorption and scattering — their combination is called “extinction” — of starlight 
are cansed by interstellar grains of all sizes. Bnt grains of sizes comparable to the wavelength 
(A) of starlight absorb and scatter photons most effectively (i.e., 27ra/A ~ 1, where a is the 
spherical radins of the grain; see Li 2009). Becanse of this, interstellar grains have long been 
known to be “snbmicron-sized”, with a canonical size of ~ 0.1 pm (i.e., a ~ Ay/27r ~ 0.1 pm 
for the visnal band A^/ = 5500 A) since their hrst detection via extinction and reddening in 
the visible (Trnmpler 1930). However, it is now well recognized that they actually span a 
range of sizes from subnanometers (i.e., angstroms) and nanometers to submicrometers and 
micrometers. 


The sub-pm-sized grain population is well constrained by the wavelength-dependent 
extinction from the near infrared (IR) at A < 3 pm to the far ultraviolet (UV) at A > 0.1 pm. 
This population, ranging from a couple hundredth micrometers to several submicrometers, 
is often referred to as “sub-pm-sized” grains or “classical” grains. 


The angstrom- and nano-sized grain population cannot be constrained even by the far- 
UV extinction. These grains (of sizes from a < 50A down to a few angstroms) are often 
referred to as “very large molecules”, “very small grains” (VSG), “ultrasmall grains”, or 
“nanoparticles”. They are in the Rayleigh regime (i.e., 2na/X -C 1) in the far-UV. On a 
per unit volume iV) basis, their extinction cross section Oext(a, A) of spherical radii a at 
wavelength A is independent of grain size, i.e., Oext(a, A)/R = Oext(A)/R. The extinction 
resulting from these grains is not sensitive to the size distribution dn/da but their total 
volume: 
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where uh (-A^h) is the hydrogen volume (column) density, and R/sg/H is the total volume 
per H nuclei of these grains. Eq.[T] demonstrates that the observed far-UV extinction is not 
able to constrain the size distribution of this grain population. On the other hand, with 
a heat content smaller or comparable to the stellar photons that heat them, these grains 
are stochastically heated by single photons and do not attain an equilibrium temperature 
(Draine & Li 2001). Upon absorption of an energetic photon of energy hv, a ultrasmall 
grain will be heated to a maximum temperature Tp determined by its specihc heat C{T) 
and hi/\ ^ C{T) dT = hu, and then rapidly radiates away most of the absorbed energy at 
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temperature Tp. Therefore, its IR emission contains crucial information about its size since 
C(T) oc (see Li 2004). Indeed, as demonstrated by Li & Draine (2001; LDOl) and Draine & 
Li (2007; DL07), the size distribution of the ultrasmall grain population is determined by the 
near- and mid-IR emission features at 3.3, 6.2, 7.7, 8.6 and 11.3/rm ascribed to polycyclic 
aromatic hydrocarbon (PAH), as well as the near- and mid-IR broadband photometry of 
IRAS and C'ORR'/DIRBE, particularly those at A < 25 /xm. 

In the UV/optical, micrometer-sized grains or “very large grains” (VLG) are in the 
geometrical-optics regime (i.e., 27 ra/A ^ 1) and their extinction is essentially wavelength- 
independent or “gray” (Li 2009). Therefore, the UV/optical extinction is not able to con¬ 
strain their quantity or size distribntion. There are several lines of direct evidence for the 
presence of very large grains in the interstellar medium (ISM). Micrometer-sized interstellar 
SiC, graphite, AI 2 O 3 , and Si 3 N 4 grains have been fonnd in primitive meteorites. They were 
identified as presolar grains of stellar origin based on their isotope anomalies (Clayton & 
Nittler 2004). Fnrthermore, the dust detectors aboard the Ulysses and Galileo spacecrafts 
have detected large interstellar grains with radii np to ~2.0/im flowing throngh the helio¬ 
sphere from the local interstellar clond (Griin et al. 1994, Kruger et ah 2007). Moreover, 
even larger interstellar grains (of radii ~ 20 //m) were detected as radar meteors entering 
the Earth’s atmosphere on solar-hyperbolic trajectories (Taylor et al. 1996, Baggaley 2000). 
More recently, Westphal et al. (2014) reported the detection of seven grains possibly of in¬ 
terstellar origin returned by the Stardust spacecraft (also see Sterken et al. 2015). These 
grains are mostly composed of Mg-rich silicates and three of them have radii > 1 /rm. 

In addition, Witt et al. (2001) pointed out that the presence of very large grains in 
the ISM could be inferred from X-ray halos since the X-ray scattering efficiency varies ap¬ 
proximately as while the differential cross section varies as a®. They examined the X-ray 
halo around Nova Cygni 1992 and found that, in order to explain the observed profile and 
intensity, the grain size distribution needs to extend to and possibly beyond ~ 2 /rm. Fur¬ 
thermore, surprisingly high-albedo scattering at near-IR wavelengths has been reported for 
several dense regions (Witt et al. 1994, Block et al. 1994, Lehtinen & Mattila 1996), sug¬ 
gesting the possible presence of a population of grains at least ~0.5 /rm in radii. Moreover, 
based on the unaccounted for 0/H abundance (Jenkins 2009, Whittet 2010), Jenkins (2009), 
Poteet et al. (2015), and Wang, Li & Jiang (2015) argued that the missing reservoir of O 
could reside on /rm-sized H 2 O grains. It is worth noting that four decades ago Greenberg 
(1974) had already pointed out that there could exist very large grains comprised of O, G, 
and N in the interstellar space. Socrates & Draine (2009) suggested that “pebble”-sized 
grains of a~ 1 mm may be detected through optical scattered light halos. Finally, we note 
that /xm-sized grains were detected in dense cloud cores through the scattering of the inter¬ 
stellar radiation at the Spitzer/IRAC 3.6 and 4.5/xm bands (i.e., “coreshine”; Pagani et al. 
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2010, Steinacker et al. 2014) @ 

In this work, we aim at constraining the quantity and size distribution of the /rm-sized 
dust population. This is achieved by htting the observed mid-IR extinction at ~ 3-8 pm (§2]) 
in terms of the silicate-graphite-PAH model together with an extra population of pm-sized 
grains (^. We present the results in §4] and discuss their astrophysical implications in ^ 
Our principal conclusions are summarized in ^ 


2. Mid-IR Extinction 

The UV/optical interstellar extinction can be characterized by a single parameter Ry 
(Cardelli et al. 1989, CCM)Jf| The UV/optical extinction can be closely htted in terms of the 
classical silicate-graphite model (Mathis et al. 1977, Draine & Lee 1984 [DL84]). In the IR at 
1 pm < A < 7 pm, this model predicts a power-law extinction curve of Ax oc (Draine 

1989). As elaborated in Wang, Li & Jiang (2014), this is too steep to be consistent with the 
subsequent observations made by the Infrared Space Observatory (ISO) and the Spitzer Space 
Telescope. Numerous observations suggest that the mid-IR extinction at 3 pm < A < 8 pm is 
flat or “gray” for both diffuse and dense environments (see Figured]), including the Galactic 
center (Lutz 1999, Nishiyama et al. 2009), the Galactic plane (Indebetouw et al. 2005, Jiang 
et al. 2006, Gao et al. 2009), the Goalsack nebula (Wang et al. 2013), and nearby star-forming 
regions (Flaherty et al. 2007). All these observations appear to suggest an “universally” flat 
mid-IR extinction law, with little dependence on environments. 


3. Dust Model 

We aim at reproducing the observed extinction from the UV/optical to the near- and 
mid-IR. We assume a mixture of amorphous silicate dust and carbonaceous dust, taking the 
size distribution functional form of Weingartner & Draine (2001; WDOl). We assume that 
the latter extends from grains with graphitic properties at radii a > 50 A, down to grains with 
PAH-like properties at very small sizes (LDOl). The WDOl model employs two log-normal 


^The presence of this type of grains in dense environments is typically interpreted as due to grain growth 
which also occurs in proto-planetary disks (e.g., see Bouwman et al. 2001, Kessler-Silacci et al. 2006, Ricci 
et al. 2010). In contrast, those p,m-sized grains discussed in this work are for the diffuse ISM as well, not 
confined to dense clouds. 

^Rv = Ay j {Ab — Ay) is the total-to-selective extinction ratio, where Ab is the blue band extinction 
at Ab = 4400 A. For the Galactic average, Ry ~ 3.1. 
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Fig. 1.— Fitting the i?v = 3.1 extinction curve from the UV/optical to the near- and mid-IR with 
(1) amorphous silicate (green short-dashed line), (2) graphite and PAHs (blue dashed line), and (3) 
sized graphite (magenta dash-dot-dotted line) with a log-normal size distribution characterized by oq ~ 
1.2/im, a ~ 0.3, and &vlg ~ 137ppm (see eq.[2]). The thick red solid line plots the model-fit which is 
the sum of silicate, graphite/PAHs, and /xm-sized graphite. The black solid line plots the sum of silicate 
and graphite/PAHs. The symbols plot the observed extinction: cyan open triangles plot the Ry = 3.1 
UV/optical/near-lR extinction, while other symbols plot the mid-lR extinction (see text). 


size distributions for two populations of PAHs which respectively peak at ao,i, ao ,2 and have 
a width of cxi, a 2 , consuming a C abundance of bc,i, &c ,2 (per H nuclei). Following DL07, 
we adopt ao,i = 3.5 A, ui = 0.40, 6c,i = 45 ppm, ao ,2 = 20 A, a 2 = 0.55, and 6c, 2 = 15 ppm. 
These parameters were constrained by the observed near- and mid-lR emission. To account 
for the mid-lR extinction, we invoke an extra population of very large grains for which we 
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Fig. 2.— Grain size distributions for (1) silicate (green short-dashed line), (2) graphite and PAHs (blue 
dashed line), and (3) /xm-sized “very large” graphite (magenta dash-dot-dotted line). 

also adopt a log-normal size distribution of peak size oq and width a: 

I dn 3 exp (—4.5cr^) fevLch^n 1 I 1 

uh da (27r)3/2 pa^a 2 a | 2 

where mp is the atomic H mass, p and p are respectively the mass density and the molecular 
weight of the dust species (p ~ 3.5gcm“^ and p ~ 172 for silicate, p ~ 2.24 gcm“^ and 
p ^ 12 for graphite), and 6 vlg is the abundance per H nuclei locked up in pm-sized grains. 
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We consider 120 wavelengths, eqnally spaced in In A, to model the extinction between 
0.125 /rm and 8 /rm. For the “observed” extinction in the wavelength range of 0.125 /im < A < 
3 /im, we take the Galactic average of Ry = 3.1 as parameterized by CCM. For the ~ 3-8 /xm 
mid-IR extinction, we hrst obtain a weighted “average” from the observed extinction shown 
in Figure [H with twice as much weight given to the diffuse sightlines toward the Galactic 
center. We then interpolate the “average” mid-IR extinction into 25 logarithmically equally- 
spaced wavelengths. 


We take the optical constants of astronomical silicate and graphite of DL84. The PAH 
absorption cross sections are taken from DL07. Following WDOl, we use the Levenberg- 
Marquardt method to minimize = Xi + Xv; where yf gives the error in the extinction 
£t: 

[In Aobs(Ai) - InAmod(Ai)]^ 


xl = ^ 


(Tf 


( 3 ) 


where Aobs(Ai) and Amod(Ai) are respectively the observed and model-computed extinction 
at wavelength A*. Following WDOl, we take the weights = 1 for 0.125 /rm < A < 
0.9/im and = 1/3 for 0.9 /rm < A < 8 /im. The “penalty” term Xv prevents the 
model-consumed G and Si abundances ([C/Hj^^gj., [Si/H]j^gJ from grossly exceeding their 
interstellar abundances: 


xl = 0.4[max(C, 1) - + 0.4[max(Si, 1) - l]^'® , (4) 

where G = [G/H]j^gj./233 ppm and Si = [Si/H]^^g^./36.3 ppm (ppm is an abbreviation of 
parts per million). 


4. Results 

In htting the far-UV to mid-IR extinction represented by 120 “data points” (see [|3]), 
we have ten parameters from the WDOl size distribution function (Gg, at,g, ac,g, Ug, (3g for 
graphite, Cg, at,s, cic,s, ccg, Ps for silicate) and three parameters from the /im-sized grain 
component. Following WDOl, we fix Oc^s = 0.1 /im0 We consider a graphite composition 
for the /im-sized grain component]^ As shown in Figured], with oq ~ 1 .2/rm, a ~ 0.3, and 
^VLG ~ 137 ppm, the observed extinction from the far-UV to the mid-IR is closely reproduced 


"^WDOl found that the fitting error only mildly varies with ac,s provided ac,s ^ 0.1 fim while it abruptly 
increases with ac,s when Oc.s >0.1 fiia. 

^In the following, unless otherwise stated, we refer “graphite/PAHs” to nano- and sub-/im-sized graphitic 
grains (including PAHs). 




(with ~ 0.23). The corresponding parameters for the WDOl size distribution function 
are: Cg Ri 5.75 x 10“^^, ag ^ —1.40, f3g Ri 0.0291, at,g Ri 0.00818 pm, and ac,g ~ 0.173 pm for 
graphite, Cg ~ 7.56 x lO”^'^, ttg Ri —2.19, /Sg Ri —0.586, and at,s ~ 0.204 pm for silicate. The 
derived size distributions are illustrated in Figure [2J 

The model requires a total silicate mass (relative to H) of M^w/M^ Ri 6.95 x 10“^, 
and a total carbonaceous dust mass (relative to H) of Mcarb/^n ~ 4.34 x 10“^, indi¬ 
cating a total gas-to-dust mass ratio of Mgag/Mdust ~ 124, where Mgas ~ Mne + Mh Ri 
1.4 Mh for He/H Ri 0.1. The required elemental depletions are [Si/Hj^^^^ Ri 40.4 ppm and 
[C/H]j^g^ Ri 362 ppm. The latter includes r: 137 ppm in pm-sized graphite grains 

and [C/H]p^jj = 60 ppm in PAHs, which respectively account for ~37.8% and ~16.6% of 
the total mass of the carbonaceous dust component, and ~ 14.6% and ~6.38% of the total 
dust mass. 

We have also tried pm-sized silicate grains. Compared with graphite, silicates are much 
more transparent at A < 8 pm. To account for the observed mid-IR extinction, one requires 
[Si/H]ypjj R 2.95 X 10^ ppm to be locked up in pm-sized silicates, far exceeding the available 
amount of [Si/H]jgj^ r 32-41 ppm in the ISM (see §5]). Also, pm-sized silicate grains exhibit 
a Si-0 resonance dip at ~ 8 pm which is not seen in the observed extinction curve. We have 
also considered pm-sized iron grains. Unless they are extremely-elongated like needles (Dwek 
2004), the mid-IR extinction requires a depletion of [Fe/Hjyp^ r 100ppm, far exceeding 
the available amount of Fe/H r 28-35 ppm in the ISM. More recently, Kohler et ah (2014) 
found that amorphous silicates with iron nanoparticles embedded could enhance the mid-IR 
extinction. Kohler et ah (2015) further found that dust aggregates could also cause increased 
mid-IR extinction. It would be interesting to see how their model extinction per H column 
compares with the observed extinction. 


5. Discussion 
5.1. IR Emission 

As shown in Figure [T|, pm-sized graphite grains are “gray” in the UV, optical, and 
near-IR. They do not absorb much in the UV/optical, and therefore by implication do not 
emit much radiation in the IR. We have calculated the temperature probability distribution 
functions of PAHs, small graphite and silicate grains of radii smaller than 250 A which 
are heated by the Mathis, Mezger, & Panagia (1983, MMP83) interstellar radiation held 
(ISRF). We have also calculated the equilibrium temperatures of silicate and graphite grains 
of radii larger than rs-/ 250 A (see Figure [3]). The resulting IR emission is shown in Figure H] 
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Fig. 3.— Equilibrium temperatures for silicate (dashed line) and graphite grains (solid line) heated by the 
MMP83 ISRF. 

and compared with that of the diffuse ISM observed by COUF/DIRBE, COBE/FIRAS 
and P/ancA:@ Our model yields a total IR intensity of ~4.55 x 10“^^ergs“^ The 

fractional contributions of silicate, graphite/PAHs, and /im-sized graphite are approximately 
27.5%, 70.0%, and 2.5%, respectively. The pm-sized graphite component is cold (T 
13.1,12.5,11.2 K for a = 1,1.2,2 pm compared to T 20K for a = 0.1pm, see Figure [3]) 
and its contribution to the overall emission is only noticeable at A > 550 pm, becomes 
signihcant at A > 1382 pm, and dominates the emission at A > 2098 pm. 

Figure 0] shows that our model fits the observed emission very well from the near-IR 
to submm: up to A < 850 pm, it is in close agreement with the COBE/FIRAS spectropho¬ 
tometry. At A > 850 pm, the model emission slightly exceeds the COBE/FIRAS data. We 
note that the Planck photometric data at A > 850 pm are also slightly higher than the 
CORF/FIRAS spectrophotometry. Our model is in excellent agreement with the Planck 
data up to A > 1382 pm, and does not overpredict until A > 2098 pm. At A = 3000 pm, 
the model emission is stronger than the COBE/FIRAS photometry by a factor of ~ 1.5 and 
exceeds the Planck data by ~ 72%. However, the far-IR dielectric functions of graphite and 


®Many Herschel observations have deen dedicated to the investigation of the dust properties in the ISM 
(e.g., see Abergel et al. 2010, Gordon et al. 2010, Molinari et al. 2010, Juvela et al. 2011, Juvela 2015). In 
this work, these Herschel data are not included since the Herschel/PACS bands are essentially covered by 
GOSE/DIRBE, while the Herschel/SPIRE bands are essentially covered by GORE/FIRAS and Planck (see 
Table 4 of LDOl and Table 6 of DL07). 




(erg s'^ sr'^ H'^) 
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Fig. 4. — Comparison of the model to the observed emission from the diffuse ISM. Green short-dashed 
line shows emission from amorphous silicate; blue dashed line shows emission from graphite/PAHs; magenta 
dash-dot-dotted line shows emission from /im-sized graphite grains. Red triangles show the model spectrum 
(red solid line) convolved with the DIRBE filters. Observational data are from DIRBE (black diamonds; 
Arendt et al. 1998), EIRAS (black squares; Einkbeiner et al. 1999), and Planck (orange circles; Planck 
Collaboration XVII 2014). 

silicate adopted here may not be known with a high precision (e.g., at A = 1000 pm the 
DL84 silicate has an opacity of Kabs ~ 0.33cm^g“^, while the silicate material measured 
by Agladze et al. (1996) had Kabs ~ 1.25cm^g“^). We also note that the Planck data at 
A = 3000 pm exceeds that of COBE/FIRA'S by ~45%. 
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Finkbeiner, Davis & Schlegel (1999; FDS) approximated the COfiF/FIRAS data in 
terms of a two-component model consisting of a warm component of temperature ^ 
16.2 K and a cold component of temperature Tc ~ 9.4 K: Ix ~ A“^-™i?A(Fvi/)+0.47 Bx{Tc). 
They ascribed the warm component to carbonaceous grains and the cold component to amor¬ 
phous silicates. We argue that while the FDS two-component model provides an excellent 
representation of the observed far-IR emission, their ascription may not be physical. As 
shown in Figure [3l to attain an equilibrium temperature of Tq = 9.4 K, silicate grains need 
to be very large, with radii a ~ 8.9 /rm. Such very large silicate grains are gray in the 
UV/optical and contribute little to the UV/optical extinction. As a result, the observed 
UV/optical extinction would exclusively rely on carbonaceous dust. However, it is known 
that neither silicates nor carbonaceous dust alone can account for the observed UV/optical 
extinction (see Li 2005a). It is more natural to attribute the warm component to a combi¬ 
nation of sub-/rm-sized silicate and carbonaceous grains and the cold component either to 
a population of pm-sized grains or alternatively, according to the two-level-system (TLS) 
model, to the low energy transitions associated with the disordered internal structure of 
sub-pm-sized amorphous grains (Meny et ah 2007). 

A cold emission component of T ~ 4-7 K was also noted in the COBE/YIRAS Galactic 
emission spectrum by Wright et ah (1991) and Reach et ah (1995). Rowan-Robinson (1992) 
explained the cold emission in terms of grains of a = 30 pm. But he did not specify the 
composition of this grain population. He simply assumed its absorption efficiency to be 
Qabsio,, A) = 1 for A < Ana and Qabs(a, A) = Ana/X for A > Ana. 

As shown in Figure lU our model SED flattens at A > 1382 pm. The overall dust opacity 
at A < 1382 pm has a power-index of /Sfir ~ 1.89. In contrast, the opacity at A > 1382 pm 
has a slightly smaller power-index of d mm ~ 1.67. The Planck observations of the Galactic 
diffuse ISM detected the flattening of the SED at long wavelengths, with /Sfir — /9mm ~ 0.15 
at A > 850 pm (see Planck Gollaboration XVH 2014). These observations also provided 
evidence for the far-IR opacity variations in the diffuse ISM, attributable to large amorphous 
carbon grains or the TLS model of disordered charge distribution in amorphous grains. 


5.2. The C/H Crisis 

The interstellar abundances of Si/H and G/H are unknown. They are often assumed 
to be that of solar ([Si/Hj^ 32.4 ± 2.2ppm, [G/H]q 269 ± 31 ppm, Asplund et al. 
2009), proto-Sun ([Si/Hj^ Ri 40.7 ± 1.9 ppm, [G/H]q Ri 288 ± 27ppm, Lodders 2003), early 
B stars ([Si/H]^ Ri 31.6 ± 3.6ppm, [G/H]^ Ri 214 ± 20ppm, Nieva & Przybilla 2012), or 
young F/G stars ([Si/H]^ Ri 39.9 ± 13.1 ppm, [G/H]^ Ri 358 ± 82ppm, Sofia & Meyer 
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2001). Our model consumes [Si/Hj^^g^ 40.4ppm which is consistent with the proto- 
Sun Si/H abundance. With the gas-phase abundance of [0/11]^^^^ 140 ppm (Cardelli 

et ah 1996) or [0/11]^^^^ 100 ppm (Soha et ah 2011) subtracted, the proto-Sun refer¬ 

ence standard leaves only ~ 148 ppm or ~ 188 ppm of C/H for the dust, respectively. With 
|C/H|,„ « 362 ppm, our model uses too much C compared to what would be available in 
the ISM. Even ignoring the pm-sized graphite component, the silicate-graphite-PAH model, 
with Ri 252ppm (WDOl), is already in “C crisis” (Snow & Witt 1995). With the 

inclusion of a population of pm-sized graphite grains, our model amplihes the so-called “C 
crisis”. 

The “C crisis” holds for all dust models, including that of Zubko et ah (2004; ZDA) 
which requires 244 ppm (and [Si/Hj^^^^ 36 ppm) and Jones et al. (2013; J13) 

which requires [C/H]j^gj. 233ppm (and [Si/Hj^^^g^ si 50ppm). Note that the mid-IR 
extinction has not been accounted for by ZDA or J13. The ZDA model produces much 
less extinction at A > 1pm than observed (see Figure 23.11 in Draine 2011). The J13 
model htted the IR extinction tabulated in Mathis (1990) which is much lower than the 
mid-IR extinction described here in §2J The total consumed C/H is expected to be similar 
to ours if the flat mid-IR extinction is to be accounted for by their models. This is indicated 
by the Kramers-Kronig relation of Purcell (1969) which relates the wavelength-integrated 
extinction to the dust quantity: Aa/Ar dA = 1.086 x Stt^F (Vdust/H) where TJust/H is the 

total dust volume per H nucleon, and F is a dimensionless factor which depends only upon 
the grain shape and the static dielectric constant of the grain material. Spherical silicates 
(F Ri 0.64) with the proto-Sun abundance contribute ~6.53 x 10“^® mag cm^ H“^ to the 
extinction integration. With Aa/Ar dA ~ 1.70 x 10“^®magcm^H“^ (see Figured]), 

one derives [C/Hj^^^^ Ri 365ppm for spherical graphite (F ^ 1.0). This is a lower limit 
since f^Ax/Niid\> Aa/Ar dA. For moderatly elongated grains, the C crisis still 

persists (e.g., with an elongation of ~2-3, one requires [C/H]j^j,^ si 275ppm of graphite of 
F ^ 1.25). 

A likely solution to the C crisis problem is that the stellar photospheric abundances 
of dust-forming elements may be considerably lower than that of the interstellar material 
from which young stars are formed. This could be caused by an incomplete incorporation 
of heavy elements in stars during the star formation process, and/or an underestimation of 
the degree of heavy-element settling in stellar atmospheres (see Li 2005b). Parvathi et al. 
(2012) derived the gas-phase C/H abundance of 16 Galactic interstellar sightlines. They 
found that [C/Hj^g^^ varies from one sightline to another, with ~ 1/3 of the sightlines having 
their gas-phase C/H abundance alone exceeding the proto-Sun C/H abundance (e.g., the 
sightline toward HD 206773 has a gas-phase C/H abundance of ~ 464±57ppm [see Parvathi 
et al. 2012] which is comparable to the total C/H abundance of ~ 462 ppm required by our 


13 


model (i.e., [C/H]^^g^ 362ppm [see [H] plus [C/Hjgg^ Ri 100ppm [see Sofia et al. 2O11])0 


5.3. An Artificial Gap in the Grain Size Distribntion? 

The size distribution derived for the /im-sized dust component is not a continuous 
extension of the sub-/im-sized component but has a gap in between (see Figure [2]). This 
seems inconsistent with the observed continuous size distribution of the interstellar grains 
entering the solar system (see Figure 2 of Frisch et al. 1999). However, the extinction 
resulting from such a continuous size distribution would be totally unlike what is actually 
observed (see Figures 3,4 of Draine 2009) @ As a matter of fact, a close inspection of the in 
situ dust size distribution of Frisch et al. (1999) reveals that it already shows a flattening off 
at a > 0.5 /im. Witt et al. (2001) also noted that the size distribution continuously extending 
to micrometers required to explain the X-ray halo around Nova Cygni 1992 would lead to 
an extinction curve with Ry 6.1. It would be interesting to see if a size distribution with 
a gap in between the sub-/im-sized grains and the /rm-sized grains like the one derived here 
could reproduce the observed X-ray halo and result in an extinction curve with Ry 3.1. 


5.4. The Origin of Very Large Grains 

The origin of /im-sized interstellar grains is not known. Theoretical calculations show 
that carbon dust can grow to sizes as large as a few micrometers in the radioactive envi¬ 
ronment of supernova ejecta, even if O is more abundant than C (Clayton et al. 1999). 
These calculations are also applicable to the formation of /rm-sized oxides within supernova 
gas having C > O. However, we do not know in the ISM how much dust originates from 
supernovae. If a substantial fraction of interstellar dust is from supernova condensates, then 
/im-sized grains may be prevalent in the ISM. 


^If we attribute the flat 3-8 /im mid-IR extinction to /tm-sized silicate or iron spheres, we would face a 
much more severe Si or Fe “crisis”: as discussed in 21 while the model consisting of /im-sized graphite grains 
requires ^60% more C/H than the proto-Sun C/H abundance, the model consisting of /im-sized silicate 
grains requires ~ 725 times more Si/H than the proto-Sun Si/H abundance, and the model consisting of 
/im-sized iron spheres requires ~ 3.0 times more Fe/H than the proto-Sun Fe/H abundance. 

®Draine (2009) added a population of large grains to the WDOl size distribution so that the new one 
approximately reproduces that of the interstellar grains detected by Ulysses and Galileo at m > 3 x 10“^^ g 
of Landgraf et al. (2000). With the inclusion of these large grains, Draine (2009) derived an extinction curve 
of Ry Ki 5.8. 
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We note that the size and quantity derived for the /xm-sized graphite dust does not rely 
on the silicate-graphite-PAH model. The /xm-sized component is inferred from the mid-IR 
extinction which is essentially independent of the observed properties (e.g., UV/optical/near- 
IR extinction, near-, mid-, and far-IR emission) previously used to constrain dust models. 
These models differ from each other in detailed composition and morphology (see Li 2004). 
But they all £t the UV/optical extinction and are deficient in comparison with the observed 
~ 3-8 /im mid-IR extinction except the WDOl Ry = 5.5 model. To account for the observed 
mid-IR extinction, a population of /rm-sized dust is required by all dust models. The size 
is inferred from a general consideration of light scattering theory. The quantity of this 
component can be inferred from the Kramers-Kronig relation of Purcell (1969). 

Finally, we note that, while the /rm-sized dust component derived here is carbonaceous 
in nature, the /xm-sized interstellar grains detected by the Stardust spacecraft are remarkably 
devoid of carbon (Westphal et ah 2014). This is probably due to the shock destruction of 
carbonaceous grains in the local interstellar cloud, as revealed by the overabundant CII in 
the interstellar gas surrounding the heliosphere (see Slavin & Frisch 2008). 


6. Summary 

The ~ 3-8 fam mid-IR extinction is observed to be flat or gray in various interstellar 
environments, including low-density diffuse clouds, translucent clouds, and dense clouds. 

It provides a sensitive constraint on very large, /rm-sized grains which are gray at optical 
wavelengths and whose existence can not be constrained by the far-UV to near-IR extinction. 

We model the flat ~ 3-8 /rm mid-IR extinction with /rm-sized grains together with a mixture 
of silicate and graphite grains of sizes ranging from a few angstroms to a few submicrometers. 

Our principal results are as follows: 

1. The observed interstellar extinction is closely reproduced from the far-UV to the 

mid-IR with [Si/H]d„^, Ri 40.4 ppm in silicate grains, Ri 362 ppm in car¬ 

bonaceous grains which include PAHs ([C/H]p^jj Ri 60ppm), submicron-sized graphite 
(C/H Ri 165 ppm), and /rm-sized graphite ([0/11]^^^ Ri 137ppm). The sizes of the 
/im-sized graphite component are modeled as a log-normal distribution peaking at 
ao ~ 1.2 jum and having a width of a ~ 0.3. 

2. Our model closely reproduces the IR emission of the diffuse ISM observed by COBE/DIRBE, 
COBE/FIRAS, and Planck from the near-IR up to A > 1382 java. With equilib¬ 
rium temperatures of T < 13 K, /rm-sized graphite grains dominate the emission at 

A > 2098 java and account for ~2.5% of the total IR emission. 



3. A carbon budget problem does arise with respect to our model when one compares 

the model-required C depletion of Ri 362 ppm with the solar and proto-Sun 

abundances of [C/H]QRi269, 288 ppm. 

4. Micrometer-sized silicate or iron grains, if present in the ISM, are unlikely responsible 
for the observed ~3-8 /xm extinction as they would require too much Si/H and Fe/H, 
exceeding the proto-Sun Si/H and Fe/H abundances by a factor of ~ 725 and 3.0, 
respectively. 
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